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We propose a new communication scheme that uses time-delayed chaotic systems with delay time 
modulation. In this method, the transmitter encodes a message as an additional modulation of the 
delay time and then the receiver decodes the message by tracking the delay time. We demonstrate 
our communication scheme in a system of coupled logistic maps. Also we discuss the error of the 
transferred message due to an external noise and present its correction method. 



PACS numbers: 05.45.Xt, 05.40.Pq 

The fact that a chaotic system, which intrinsically 
posses unpredictability and a broad band spectrum, can 
be synchronized 0, U Q has led people to assume that 
a chaotic system would provide better security in com- 
munication than other cryptographic schemes previously 
proposed. Since the first demonstration of its possibil- 
ity in an electronic circuit chaos communication has 
been extensively investigated and correspondingly var- 
ious methods for masking the message have been de- 
veloped It was shown, however, that the message, 
when masked by a chaotic signal from a low-dimensional 
chaotic system, can be extracted 0- Further, it was 
reported that even when the message is masked by a hy- 
perchaotic signal, it can be extracted by using nonlinear 
dynamic forecasting as far as the local dynamics does not 
reflect more complicated dynamics significantly 0. 

A cure for the weakness of the conventional chaos com- 
munication was to develop methods u sing synchroniza- 
tion of time-delayed systems [1 II E3 El El E3 De- 
spite of the small number of physical degrees of freedom, 
the time-delayed system has the advantage of possess- 
ing the property of high-dimensional hyperchaos which 
can be easily implemented electronically 9]. This made 
the communication using time-delayed systems attract 
much attention [9L llfjj. However, it was again reported 
that delay time t can be also detected by analyzing the 
transmitted signal and that in such a case the recon- 
structed phase space of the time-delayed system collapses 
into low-dimensional manifold E3- O n this discovery, it 
was demonstrated that the message masked by the sig- 
nal of time-delayed system can be extracted even in the 
presence of message signal of small amplitude EH- So 
far almost all of the methods proposed for chaos com- 
munication have been broken by their successive counter 
examples. This fact has spurred a debate on the general 
assumptions on a chaotic system. 

Meanwhile, a time-delayed system with delay time 
modulation (DTM) in which the delay time is driven by 
chaotic or stochastic signal for the purpose to make the 
delay time undetectable was introduced |16| . It was ana- 
lyzed that a phase space reconstruction is hardly possible 
in that case E3 because DTM significantly increases the 



complexity of an attractor and it renders the delay time 
undeterminable. The report that robust synchronization 
can be established between two coupled chaotic systems 
with DTM calls attention to a time-delayed system with 
DTM as an ideal candidate for communication. In the 
synchronization with DTM, the system consists of two 
parties. One is a transmitter, which is a time-delayed 
system with delay time modulation: 

x = f{x(t),x(t - r)), 

r = g(x(t),t), (1) 

where g(x(t),t) is the modulation function. The other is 
a receiver which is a time-delayed system driven by the 
signal x(t — t) from the transmitter: 

x' = f(x'(t),x(t-r)). (2) 

Even though the modulation works only at the transmit- 
ter side, it has been found, the receiver is synchronized 
with the transmitter by the signal x(t — r) |16| . 

In this Letter, we propose a communication method 
that hides a message into the delay time of time-delayed 
systems with DTM. In our method the transmitter en- 
codes the message as an additional modulation of the 
delay time by which we overcome the drawbacks in the 
previous schemes. The receiver decodes the message by 
identifying the modulated delay time with a delay buffer. 

First, we shall describe the proposed scheme concep- 
tually and demonstrate it in two coupled logistic maps. 
The transmitter is given by: 

x = f(x,x(t - f)), 
t = g(x), 

f = T + m{t). (3) 

Where we call r the bare delay time which is modulated 
by the function g{x) and f the genuine delay time which 
includes the message m(t). Here we start by supposing 
that the g(x) depends only on the state variable x, and 
that g(x) is announced in public. Even if random people 
know g(x), they can not get bare delay time r because for 
that they also need to know the synchronized state vari- 
able x. Now we consider the case that the transmitter 
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and receiver are already synchronized after some tran- 
sient time. For communication, the transmitter sends 
x* = x(t — f) which includes the encoded message. We 
emphasize here that our encoding method is fundamen- 
tally different from the previous ones. Up to now, in 
the conventional methods the encoding of the mes- 
sage was usually done by perturbing the real trajectory 
of the chaotic signal such that x = x + m. Accordingly, 
the message could be extracted by identifying the fluctu- 
ations of the reconstructed trajectory on phase space or 
return map Q. Moreover, because the message plays a 
role of effective noise which degrades the quality of syn- 
chronization, the system should have fast relaxation to 
the attractor to keep the communication quality, which 
potentially increases vulnerability pj, Ej . On the other 
hand, in our scheme, the message is encoded as a modu- 
lation of the delay time such that x — (t+tm)) which 
is just the temporal shifting of the original trajectory de- 
pending on the message to. Therefore the amplitude of a 
message need not to be small. Furthermore, since the de- 
lay time including the message changes the genuine char- 
acteristics of the attractor like the embedding dimension 
and the complexity through the modulation of the delay 
time , the message is constituted into the attractor 
itself rather than a perturbation. 

To decode the message, the receiver should preserve 
the history of its own trajectory for the time interval 
[t — r m , t] in the delay buffer which is denoted by the 
symbol {a^[i]}t_ T j where r m = maxr(i), i £ [0,N]. 
Here N — \j m /St] where St is a sampling step, and [x] 
is the largest integer less than x. For tracking the delay 
time, the receiver defines the delay identification measure 
like this: 



M(i,e) = e-\x*- {x'[i\} 



(4) 



where e is the predefined threshold for the identification 
and | • | is the absolute value. The receiver can find the 
value of f by finding the index i* which maximizes the 
delay identification measure such that M(i*, e) > M(i, e) 
for i S [0, N — 1]. Accordingly, the two delay times, i.e., 
bare and genuine delay times, can be obtained following 
the procedure: 



x> = f(x',x(t-f)), 
r' = g(x% 



N 
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(5) 



and then the message hidden in the delay time can be 
decoded at the receiver side such that: 



m\t) 



(6) 



Since we have supposed that, the two systems in our 
scheme are synchronized, the state variables and the de- 
lay times are coincided with each other, i.e., x = x' , 
t = t' and f = f'. Thus the decoded message at the 
receiver side m'(t) is equal to the original message m(t). 
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FIG. 1: Temporal behaviors of the coupled systems in Eq. 
(7) and (8), when A = 100. (a) The modulated delay time; 
The difference of two state variables: x — x' (b) below the 
threshold a = 0.28 (c) above the threshold a — 0.33. 



In conventional schemes, since the message is encoded 
as a deviation from the synchronized trajectory, the en- 
coded message can be partly decoded through the syn- 
chronized windows in such cases as a different system is 
located on the intermittently synchronized phase. How- 
ever, in our scheme, even if a different system is located 
on the intermittently synchronized phase, the message 
is completely unknown, because the message is not a 
deviation from the synchronization manifold. Another 
strength of this scheme is that the delay time is repre- 
sented by the integer even if one use a chaotic flow for 
transmitter and receiver, because delay time is actually 
an indicator of a position of the previous state. For that 
reason, digital data can be directly encoded and decoded 
in this scheme. 

For demonstration, we consider two logistic maps: 



X n+1 = Xx n (l- 

T = [Ax n ], 

f = t + to mod t„ 



Transmitter 



(7) 



where we take g(x) — [Ax] and 



Xx' n (l 
[A<], 



X n)i 



Receiver 



(8) 



to = r — r 



wherex„ = (1— a)x n +ax n -f andec^ = (1— a)x' n +ax n -f 
and we take A = 4.0. Here i S [0, r m — 1] and i* is the 
identified index which maximizes the delay identification 
measure (Eq. (4)). Figure 1 shows the temporal behav- 
iors of the two coupled systems with a null message, i.e., 



= 0. While the difference 



between the two state 



variables shows the intermittent chaotic bursting below 
the synchronization threshold (Fig. 1(b)), it converges 
fast to the synchronization manifold x — x' — above 
the threshold (Fig. 1(c)). 



(a) 




FIG. 2: (a) The conditional Lyapunov exponent A c as a func- 
tion of A and a. (b) The contour plot in the (A, a) space. 
Here the colors indicate the different values of A c . 
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FIG. 3: (a) The transmitted signal including the message m 
at a = 0.7, A = 100, and e = 10" 8 ; (b) The autocorrelation 
of the transmitted signal, which shows 8 function shape; (c) 
The probability distribution function of the transmitted signal 
x n -f ; (d) The original message at transmitter side; (e) The 
decoded message at receiver side. 



To analyze the critical behaviors of the two coupled 
chaotic systems we find the difference motion such that: 
AI n+ i = X(l-a)(l-(x n +x' n ))AX n where AX n = x n - 
x' n . The conditional Lyapunov A c , which determines the 
synchronization threshold, can be found by following the 
standard procedure: A c = lirrin-foo - log \AX n /AX \. 
Figure 2 shows conditional Lyapunov exponent as a func- 
tion of coupling strength a and modulation amplitude 
A. One sees that the conditional Lyapunov exponent be- 
comes negative above a = 0.32 for all A s [3, 100] and 
it means two coupled systems are synchronized in that 
regime. On the other hand, if the modulation amplitude 
is the relatively small, i.e., A < 3, synchronization is es- 
tablished in the stronger coupling, i.e., a > 0.48 and some 
synchronization islands appear in the regime, A £ [3, 18] 
and a < 0.32 (see the contour plot in Fig. 2 (b)). 
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FIG. 4: The decoded messages for different values of the 
threshold e in the presence of noise of max(£ n ) = 10 -6 . (a) 
The original message; (b) The decoded message with e = 
2.5 x 10~ 5 . (c) With e = 2.0 x 10~ 5 ; (d) With e = 1.5 X 10~ 5 . 
The arrows in (c) and (d) indicate the bit-flipped events in 
the decoded message 



Figure 3 shows the temporal behavior and statistical 
properties of the transmitted signal with an example of 
the message transfer. The transmitted signal x n -f is 
shown in Fig. 3 (a) and the autocorrelation function of 
the signal is presented in Fig. 3 (b). Figure 3 (c) shows 
the probability distribution function P of the transmit- 
ted signal x n —r , which is a normalized histogram of the 
projected trajectory onto the x n -f axis in n versus x n -f 
plot (Fig. 3 (a)). Since the autocorrelation function is 
^-correlated, one can see that the system looks random 
to an eavesdropper. Fig. 3 (d) is the original message en- 
coded at the transmitter and (e) is the message decoded 
at the receiver. Consideration of the effects of noise in the 
transmission channel is essential in regard to real applica- 
tion and implementation. The noise in the transmission 
channel induces the distortion at receiver side such 
that: 

x* = x* + £„, x' n = x' n + 0K„). (9) 

The distortion is propagated into the delay identification 
measure such that M(i, e) = e — \x* — {^'[«]}"_ Tm |- So 
there exists a possibility that the identified index i* can 
be determined incorrectly due to the external noise. Fig- 
ure 4 shows how the external noise has an effect on the 
message transfer in our scheme. The original message is 
presented in Fig. 4 (a) and the decoded messages with 
different values of the threshold e for delay time identifi- 
cation are presented in Fig. 4 (b)-(d). One can see that 
the transferred bits are intermittently flipped as shown 
in Fig. 4 (c) and (d) (see the arrows). To clarify the role 
of the noise in our scheme and the relationship with 
the threshold e, we evaluate the average decoding error 
which is defined by the number of flipped bits divided 
by the total number of transferred bits. Figure 5 (a) 
shows the decoding error as a function of the threshold 
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FIG. 5: The average decoding errors for different noise ampli- 
tudes as a function of the threshold e, where the points below 
10 _6 '° are zero. Each datum point is the result of 10 6 bit 
transfer. Circles: max(£) = 10~ 6 ; Squares: max(£) = 10 -7 ; 
Diamonds: max(£) = 10 -8 . (a) Without the error correction; 
(b) With the error correction. See Ref. [l8| for details of the 
error correction method. 



e in the presence of external noise. The decoding error 
is decreased as the threshold is increased. Figure 5 (b) 
shows the data on the application of the error correction 
method • With the proper threshold, the receiver can 
decode the message completely even in the presence of 
external noise. 

It is worth discussing the possibility of eavesdropping 
for the proposed scheme. The eavesdropper can accumu- 
late the transmitted signal x(t — f), which is the tem- 
porally shuffled signal of the synchronized state variable 
x(t). Accordingly the eavesdropper can not reconstruct 



the phase space correctly 01 ; because the information of 
the phase space is preserved only when the temporal or- 
dering of the signal is kept. Even if one succeeds in recon- 
structing the phase space, the message is not separable 
as mentioned above. For the eavesdropper to identify the 
genuine delay time f , the construction of the delay buffer 
is the most essential procedure. The eavesdropper, how- 
ever, always gets the fake delay buffer {i[i]}(_ Tm in which 
the temporal order is mixed. Accordingly the eavesdrop- 
per read the incorrect genuine delay time f from it. Fur- 
thermore, since the synchronized state variable x is never 
exposed outside the transmitter or the receiver in our 
scheme, the eavesdropper can not find the bare delay 
time r also, which is the function of synchronized state 
variable x. 

In conclusion, we have proposed a new communication 
scheme which enables one to encode and decode the de- 
sired message using the time-delayed system in the pres- 
ence of delay time modulation. Our scheme is fundamen- 
tally different from the conventional ones in that, in our 
scheme, since the message is encoded as a modulation of 
the delay time, it does not degrade the quality of the syn- 
chronization; in the conventional ones the message is en- 
coded as a perturbation of the real trajectory of a chaotic 
system and so the quality of synchronization is degraded, 
which eventually increases the potential vulnerability of 
whole communication systems. Also the message is de- 
coded just by comparing the transmitted signal with the 
history of the trajectory stored in the delay buffer. We 
have also shown that the scheme works even in the pres- 
ence of external noise just with a simple algorithm for 
the error correction. We expect our scheme can be used 
to implement the real communication system with better 
performance and enhanced security. 



Acknowledgments 

This work is supported by Creative Research Initia- 
tives of the Korean Ministry of Science and Technology. 



[1] H. Fujisaka and T. Yamada, Prog. Theor. Phys. 69, 32 
(1983); V.S. Afraimovich, N.N. Verichev, and M.I. Rabi- 
novich, Radiophys. Quantum Electron. 29, 747 (1986); 
L.M. Pecora and T.L. Carroll, Phys. Rev. Lett. 64, 821 
(1990). 

[2] S. Boccaletti, J. Kurth, G. Osipov, D.L. Valladares and 
C. Zhou, Physics Reports 366, 1 (2002). 

[3] A. Pikovsky, M. Rosenblum, and J. Kurths, Synchro- 
nization A universal concept in nonlinear science, CAM- 
BRIDGE UNIVERSITY PRESS, 2001. 

[4] K.M Cuomo and A.V. Oppenheim, Phys. Rev. Lett. 71, 
65 (1993). 

[5] L. Kocarev and U. Parlitz, Phys. Rev. Lett. 74, 5028 
(1995); J. H. Peng, E. J. Ding, M. Ding, and W. Yang, 
Phys. Rev. Lett. 76, 904 (1996); J. H. Xiao, G. Hu, and 



Z. Qu, Phys. Rev. Lett. 77, 2818 (1996); S. Boccaletti, A. 
Farini, and F. T. Arecchi, Phys. Rev. E 55, 4979 (1997); 
K. Murali and M. Lakshmanan, Phys. Rev. E 56, 251- 
255 (1997); Z. Liu, S. Chen, and B. Hu, Phys. Rev. E 
59, 2817 (1999); S. Sundar and A. A. Minai, Phys. Rev. 
Lett. 85, 5456 (2000); C.-M. Kim, S. Rim, and W.-H. 
Kye, Phys. Rev. Lett. 88, 014103 (2002); S. Wang, J. 
Kuang, J. Li, Y. Luo, H. Lu, and G. Hu, Phys. Rev. E 
66, 065202 (2002); R. Mislovaty, E. Klein, I. Kanter, and 
W. Kinzel, Phys. Rev. Lett. 91, 118701 (2003). 
[6] G. Perez and H.A. Cerdeira, Phys. Rev. Lett. 74, 1970 
(1995). 

[7] K.M. Short and A.T. Parker, Phys. Rev. E 58, 1159 
(1998); D.-U. Hwang, C.-M. Kim, and Y.-J. Park, J. Ko- 
rean Phys. Soc. 42, 8 (2003). 



■5 



[8] K. Pyragas, Phys. Rev. E 58, 3067 (1998); R. He and 
P.G. Vaidya, Phys. Rev. E 59, 4048 (1999). 

[9] B. Mensour and A. Longtin, Phys. Lett. A 244, 59 
(1998). 

[10] J. -P. Goedgebuer, L. Larger, and Henri Porte, Phys. Rev. 

Lett. 80, 2249 (1998); L. Yaowen, G. Gguangming, Z. 

Hong, and W. Yinghai, Phys. Rev. E 62, 7898 (2000); 

V.S. Udaltsov, J.-P. Goedgebuer, L. Larger, and W.T. 

Rhodes, Phys. Rev. Lett. 86, 1892 (2001). 
[11] J.K. Hale, Theory of Functional Differential Equations, 

Springer- Verlag, Berlin, (1997) and references therein. 
[12] T. Heil, I. Fischer, W. Elsasser, J. Mulet, and C. R. Mi- 

rasso, Phys. Rev. Lett. 86, 2001 (795). 
[13] J. Fort and V. Mendez, Phys. Rev. Lett. 89, 178101 

(2002). 

[14] R. Hegger, M.J. Biinner, H. Kantz, and A. Giaquinta, 



Phys. Rev. Lett. 81, 558 (1998). 
[15] C. Zhou and C.-H. Lai, Phys. Rev. E 60, 320 (1999); 
B.P. Bezruchko, A.S. Karavaev, V.I. Ponomarenko, and 
M.D. Prokhorov, Phys. Rev. E 64, 056216 (2001); V.I. 
Ponomarenko and M.D. Prokhorov, Phys. Rev. E 66, 
026215 (2002). 

[16] W.-H. Kye, M. Choi, M.-W. Kim, S.-Y. Lee, S. Rim, C- 
M. Kim, and Y.-J. Park, Phys. Lett. A 322, 338 (2004). 

[17] W.-H. Kye, M. Choi, S. Rim, M.S. Kurdoglyan, C.-M. 
Kim, and Y.-J. Park, Phys. Rev. E 69, 055202(R) (2004). 

[18] As a simplest implementation of error correction, we have 
introduced a redundance of the transferred data. At the 
transmitter, a bit is represented by five duplicated bits, 
e.g., "1" is represented by "11111". At the receiver, the 
bit is read as the value that is duplicated more than two 
times. 



